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ABSTRACT 


The  hydrodynamic  equations  of  a  hot  non-uniform 
plasma  are  solved  numerically  in  slab  geometry  to  obtain  the 
resonant  frequencies  and  associated  wave  functions.  The 
splitting  of  the  various  resonances  is  Bhoun  to  depend  on  the 
parameter:  (slab  thickness  «  Debye  length). 


Resonance  Oscillations 


In  a  Hot  Kon-Uniform  Plasma 
Jerald  V.  Parker 


Until  recently  the  large  splittings  of  the  Dattner  resonances  have 
remained  unexplained,  although  a  small  splitting  can  be  shown  to  result 
from  non-zero  temperatures'1'. 
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In  a  recent  letter  Weissglas  has  shown  that  results  in  qualitative 
agreement  with  observations  can  be  obtained  by  assuming  a  non-uniform 
electron  density.  In  order  to  facilitate  comparison  of  the  theory  with 
experiment  It  is  necessary  to  calculate  the  frequencies  of  the  resonances 
for  some  approximately  correct  electron  density  distributions. 

The  equations  used  below  to  describe  the  behavior  of  a  non-uniform 
plasma  in  the  collisionless  approximation  are  the  linearized  hydrodynamic 
equations  and  Maxwell's  equations.  Since  experiments  have  generally  been 
conducted  In  plasmas  for  which  one  free  space  wavelength  is  large  com¬ 
pared  to  the  dimensions  of  the  plasma,  one  can  make  a  quasi- static 
analysis .  > 

The  equations  are  then 
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where  n f(r)  is  the  steady  state  electron  density  and  where  it  has  been 


2. 


assumed  that 

n  -  nQf(r)  +  ^  t (0)  -  1 
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Since  these  waves  will  propagate  adiabatieally  we  can  relate  and 
as  follows: 

Px  =  Tin  Dj.  (4) 

Combining  equations  (1)  through  (4),  assuming  an  e’1*0*'  dependence,  and 
letting  t1  •  -  VjZ^ 


are  the  plasma  frequency  and  the  Debye  length  at  the  point  .r  =  0  . 
Despite  the  complexity  of  this  equation  the  calculations  progress  quite 
easily  if  one  chooses  to  york  in  plane  geometry. 

If  one  considers  the  one-dimensional  problem  illustrated  in  Figure  1, 
the  following  simplifications  can  be  made.  First  only  exists  which 
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Figure  1. 


reduces  the  order  of  the  equation  by  one;  second,  the  resulting  equation 
can  be  Integrated  once  to  yield 
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Using  the  boundary  condition  that  normal  vanishes  at  the  edge  of  the 

plataa,  one  can  evaluate  the  constant  K  . 


K  =  -  ——x  E,  (d)  (7) 

The  final  slapllfi  cation  results  If  one  changes  to  the  dimensionless 
variable  s  =  z/d  and  recognizes  that  a  resonance  in  this  plane  geometry 
is  characterized  by  the  vanishing  of  the  external  field  so  that  one  need 
solve  only  the  hoaogeneous  equation 
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The  density  functions  f(z)  used  in  these  calculations  are  taken 
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from  the  work  of  S.  Self  and  represent  theoretical  curves  for  a  low 
density  plasma.  They  are  shown  in  Figure  2.  The  small  slope  discon¬ 
tinuities  are  due  to  the  approximation  used  in  the  calculations. 

The  results  of  numerical  solution  of  this  equation  for  several 
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values  of  the  parameter  (r— )  are  shown  in  Figure  3.  The  circle 

u 

points  are  the  results  of  numerical  solution  of  the  differential  equation. 
The  triangle  points  were  found  using  the  WKB  method  since  the  differen¬ 
tial  equation  is  difficult  to  solve  for  very  short  wavelengths.  The 
dotted  lines  connecting  the  points  are  included  to  show  which  points 
belong  to  the  same  mode  and  to  indicate  the  general  trend  of  the  fre¬ 
quencies  but  do  not  represent  actual  data.  The  solid  line  shows  the 
ratio  of  mean  density  to  the  square  of  the  frequency  of  the  lowest  mode 
n  f  /  os2 

In  general  it  may  be  seen  that  the  splitting  of  the  resonant  fre¬ 
quencies  is  dependent  on  the  parameter  (d/x^)  and  the  splitting  is 
largest  for  small  values  of  this  parameter. 

Figures  4  and  5  show  how  this  strong  splitting  is  induced  by  the 
non-uniform  plasma  density.  Both  figures  show  the  electric  field  at 
resonance  for  various  different  modeB.  Figure  4  shows  the  modes  o>n> 
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for  a  fixed  value  of  (d/x^)  .  At  a  given  frequency  00 
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that  portion  of  the  plasma  for  which  cd  >  — —  (z  >  z  ) 
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propagate  plasma  waves,  while  the  remainder  cannot.  As  co  is  increased 


the  propagating  region  becomes  longer  and  successively  more  half 


5 


wavelengths  of  the  v*va  oan  be  fit  Into  it. 

&  2 

In  Pigure  5  the  mode  ui  is  shown  for  various  values  of  (~)  . 
d  2 

One  can  see  that  as  (— )  Increases  the  wavelength  of  the  oscillation 

u 

decreases.  Since  the  thickness  of  the  sheath  and  the  wavelength  are 
proportional  to  X^  the  resonant  frequency  of  the  mode  u)q  ,  which  can 
propagate  primarily  in  the  sheath  region,  tends  to  remain  constant. 

It  is  conceivable  that  this  effect  might  cause  one  or  more  reson¬ 
ances  to  remain  at  frequencies  below  the  lowest  plasma  frequency 
n  f(d)  e2  /  2 

-p-  — -  (see  Fig.  2  for  (d/X^)  -  a)  in  the  limit  X^  -*  0  .  This 

o 

would  mean  that  the  resonances  would  remain  split  even  in  the  limit  of 
zero  temperature  and  that  the  curves  in  Figure  3  would  tend  to  a  limit 

>  .1  as  (^-)  -►  oo  .  This  does  not,  however,  seem  to  be  the  case 
l) 

although  one  cannot  be  certain  of  extrapolations  based  on  ouch  limited 
data. 

I  wish  to  express  my  appreciation  to  R.  W.  Gould  and  to  the  Office 
of  Naval  Research  and  the  National  Science  Foundation  for  supporting  this 
research. 
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